Critical exponents of directed percolation measured in spatiotemporal intermittency 
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A new experimental system showing a transition to spatiotemporal intermittency is presented. 
It consists of a ring of hundred oscillating ferrofluidic spikes. Four of five of the measured critical 
exponents of the system agree with those obtained from a theoretical model of directed percolation. 
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I. INTRODUCTION 

Intermittency was observed in hydrodynamics as a pre- 
cursor to turbulence (see e.g. Q). For dissipative dynam- 
ical systems without spatial degrees of freedom intermit- 
tency was first modeled by ordinary differential equations 
and iterative maps (2). Spatiotemporal intermittency 
(STI) is a further development of this concept for spa- 
tially extended systems and was introduced by Kaneko 
H in the mid 80's Q]. It is characterized by patches 
of ordered and disordered states fluctuating stochasti- 
cally in space and time. This behaviour has been ob- 
served in many systems. Theoretical approaches have 
been made in a large variety of systems ranging from par- 
tial differential equations, such as the damped Kuramoto- 
Sivashinsky equation fi, |(J or the complex Ginzburg- 
Landau equation @J8T^|, over stochastical partial dif- 
ferential equations (l(J to coupled map lattices (CML) 
H |ll| and probabilistic cellular automata Q. 

In 1986 Pomeau [|l3| suggested that the onset of chaos 
via STI might be analogous to directed percolation (DP) 
p| . Such processes are modeled as a probabilistic cellu- 
lar automaton with two states per site, associated with 
the laminar and chaotic patches in the case of STI. One 
of the main features of DP is the presence of an absorb- 
ing state, which corresponds to the laminar state. The 
absorbing state prevents the nucleation of chaotic do- 
mains within laminar domains. DP model predicts some 
universal properties of STI. In particular, the fraction of 
chaotic domains is expected to grow with a power law 
e^, where e measures the distance from threshold. The 
correlation length decreases with e~ Vs , as does the corre- 
lation time with e _zyt , and the critical distribution of the 
laminar lengths is determined by / _/is and by t _/it for the 
laminar times. 

Some experimental tests of the conjecture |fL3[ in quasi 
one-dimensional systems have been made [f7[ |15|, |l6| , |l7| , 
|l8| , |l9| , p0| , |2l| , |22]] . A short summary of these and the 
relevant exponents |23|, Q for the comparison with the 
DP-model are given in Table |. However the statement 
that "...there is still no experiment where the critical be- 
haviour of DP was seen" by Grassberger |25| still seems 
to be true. Thus in this paper a new experimental ap- 
proach to this old problem is presented. 

The system introduced consists of a ring of ferrofluidic 
spikes excited by an external magnetic field. Ferrofluids, 



also called magnetic fluids, are a collodial suspension of 
ferromagnetic nanoparticels. The fluid is superparamag- 
netic Q . The idea for using this fluid was motivated by 
the fact that a single peak of ferrofluid can show chaotic 
oscillations under external driving of a magnetic field 
p7| , p8[ . We introduce here a system where about 100 
of these oscillating peaks are coupled by magnetic and 
hydrodynamic interaction. They exhibit changes in peak 
height of about 10 % and variations in wavelength A of 
about 50%. This system is advantageous because of its 
short response times and the easy control of the excita- 
tion. 

The article is organized as follows: In Sec.|| we de- 
scribe the experimental setup, the procedure of the mea- 
surement and the data extraction methods. In Sec.[n| 
the quantitative results are presented, i.e. the critical 
exponents (3, v Sl v t , p s and Finally in Sec.[[y| the 
results are discussed and an outlook to further investiga- 
tions is given. 



II. EXPERIMENT 

The experiment is based on the Rosensweig instabil- 
ity [|6| . This instability is observed in a horizontal layer 
of magnetic fluid, when a threshold of the vertically ori- 
ented magnetic field is surpassed. The flat surface be- 
comes unstable and a pattern of liquid spikes emerges. 
In case of an inhomogenous field the wavelength of this 
pattern can be controlled by the magnetic field H. The 
wavelength scales with the gradient of the magnetic field 
divided by a basic field at the undisturbed surface of the 
fluid. A larger gradient emphasizes small peaks and sup- 
presses large ones. Thus, inorder to generate as many 
spikes as possible, the setup of the experiment consists 
of a cylindrical electromagnet with a sharp edge (Fig.J^). 
The magnetic fluid is trapped by the inhomogenous mag- 
netic field at this edge of the magnetically soft iron core. 
In that way the 40 mm diameter of the pole shoe supports 
a ring of up to 130 spikes of magnetic fluid as indicated 
by the picture in Fig.|l|. 

The magnetic fluid used in this experiment is EMG901 
from Ferrofluidics™, a fluid based on magnetite 
Fe 3 4 with isoparafin as carrier fluid. At 20° C EMG901 
has a density p = 1.53 kg • m -3 , a surface tension a = 
29 mN • m _1 , a dynamical viscosity rj = 25mPas and a 
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TABLE I: Experiments and results in quasi onedimensional. 
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FIG. 1: Sketch of the experimental setup. 



susceptibility £ = 3. 

The ring of spikes is recorded with a CCD-camera 
mounted above the pole shoe. The electromagnet con- 
sists of a bias coil and an excitation coil. The bias coil is 
provided with a direct current of / = 1.0 A to keep the 
magnetic fluid in its place. The excitation coil is driven 
by an alternating current, phase- locked with the camera 
frequency, providing a stroboscopic jitter free recording 
on long timescales. The alternating current can be ad- 
justed between and 4.1 A. In this interval the number 
of spikes ranges from 60 to 130. For the amplitudes used 
in the experiment 108 spikes have been observed. 

To keep the viscosity and the surface tension constant 
the fluid is temperature controlled to 12.5 ± 0.03 °C by 
cooling the pole shoe. To prevent the evaporation of the 
isoparafin, the volume around the egde of the pole shoe 
is sealed with a glass plate to provide long term stability. 

The spatiotemporal behaviour is investigated by ob- 
serving the ring with the CCD-camera. To extract the 
wavelength and amplitude of the spikes as a function of 
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FIG. 2: Space time plots of the different states of the system, 
(a) /ex = 2.8 A: laminar state; (b) 7 e x = 3.0 A: spatiotempo- 
ral intermittency; (c) 7 e x = 3.6 A: chaotic state. 500 excita- 
tion periods are shown. The position x is normalized over the 
size of the ring. 



space, we define a ring of interest around the center of 
the pole shoe which is covering the ring of spikes. It is 
divided into 1024 segments. The average of the grey val- 
ues within each segment represents the amplitude. That 
way we get a spatial resolution of around 10 segments 
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FIG. 3: Complexification of the data (see text). The left hand side displays the data in real space ((a),(c)), the rigth hand side 
the corresponding graphs in Fourier space ((b),(d)). 



per peak. This reduction of the 2D image to a single line 
scan can be done in real time with a frequency of 12.5 Hz. 

In Fig.||(a) 500 of such scans of a laminar state are 
shown in space and time, where dark regions correspond 
to high amplitudes. The driving frequency of the exci- 
tation field is / ex = 12.5 Hz as mentioned above. The 
period r = l// ex is used to scale the time. Due to this 
stroboscopic recording the oscillations of the spikes can- 
not be seen. 

For the measurements a current of 1 A is applied to 
the bias coil and a constant volume of the magnetic fluid 
is dropped on the egde of the pole shoe. After a wait- 
ing time of 2 h thermal equilibrium is reached. Then 
the second coil is provided with a sinusodial excitation 
signal of an amplitude of 7 ex = 4.05 A, driving the sys- 
tem into the fully chaotic regime similar to the one indi- 
cated in Fig.||(c). That state is the basis for a quench 
to a lower excitation value corresponding to the STI 
regime (Fig.|](b)), which is then analyzed subsequently. 
The recording of the data starts after a waiting time of 
~ 1800t, when the transients following the quench have 
died out. The data are recorded for 2000 excitation pe- 
riods r. For higher excitation amplitudes the laminar 
state (Fig. || (a)) becomes intermittent in space and time 
(Fig. @(b)) and eventually chaotic (Fig.||(c)). 

For a quantitative analysis of the transition a reliable 
detection of the regular domains is fundamental. Differ- 
ent criteria for the distinction of the regular and chaotic 
domains have been proposed [^5j -[§2|. Because of the 
strong variations of the wavelength of our system we have 
selected a wavelength criterion, which is based on the lo- 
cal wavelenght A(x, t). To obtain A(x, t) we use a method 
called complexification. By this method a zero imaginary 
part is added to every real value (see Fig.||(a)). Then the 
data are transformed to Fourier space, where the ampli- 
tudes corresponding to negative wavenumbers are elim- 
inated (Fig.g(b),(d)). With a backward transformation 
to real space, every value has a non zero imaginary part 
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FIG. 4: (a) Local wavelength A calculated from Fig,pKc). 
(b) Relative change of local wavelength A. The dashed line 
corresponds to the threshold of 0.01 taken from Fig.|i|. 



(Fig. ||(c) dashed line). A is subsequently calculated as 
the phase difference between neighbouring values in real 
space (Fig.^(a)). In the last step the relative change of 
the local wavelength 



t + 1 



Atl 



At 



(1) 



is calculated. To get a clear distinction between reg- 
ular and irregular domains, changes in A which are 
larger than 0.01 are counted as irregular, whereas smaller 
changes belong to regular domains (Fig.^(b)). The 
threshold value A c is derived from the distribution of 
A for the fully laminar case which is presented in Fig.|| 
by a solid line. There are no larger variations than 
1%. Smaller values of A are artefacts of the recording 
technique and thus are suppressed. For comparison the 
dashed (dotted) lines give the distribution of A for the in- 
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FIG. 5: Distributions of the relative change of wavelength A 
for three different 7 e x: solid line 7 e x = 2.8 (A) correspond- 
ing to Fig.^(a), dashed line I ex = 3.0(A) corresponding 
to Fig.g(b) and dotted line 7 e x = 3.5(A) corresponding to 
Fig. Q(c). The cutoff of the distribution for the fully laminar 
state gives the threshold of the binarization A c = 0.01, as 
depicted by the arrow. 



termittent (chaotic) states respectively. The calculations 
of the exponents are robust to changes of the threshold of 
up to 50%. This variation of the binarization threshold 
corresponds to changes in the resulting critical exponents 
within the statistical errors. In Fig. |] we demonstrate the 
application of the above described procedure to the inter- 
mittent data. Fig.^(a) gives the raw data, (b) displayes 
the local wavelength and (c) the relative change A after 
binarization. The quantitative analysis of STI described 
below is based on these binarized information. To get a 
better signal-noise-ratio we average the results over six 
independent runs of the experiment, which include six 
refills of the apparatus with fresh fluid. 



III. EXPERIMENTAL RESULTS 

At low excitation amplitudes I ex <C 3.0 A the system is 
completely regular (Fig.|| (a)) showing 108 spikes. Slight 
spatial variations of the wavelength A < 0.01 remain 
constant in time. In Fie;. | (b) at /ex = 3.0 A irregular 
fluctuations are apparent, which we consider as a mani- 
festation of STI. A fairly clear distinction between regular 
and irregular domains can be made in this image even by 
naked eye. Further increase of 7 ex leads to a spreading of 
the irregular domains engulfing the regular regions, until 
finally the whole system is chaotic (Fig.|| (c)). 

As an order parameter for STI we take the time- 
averaged chaotic fraction 7, which is the ratio of chaotic 
regions to the length of the system. This ratio is averaged 
over the 2000 excitation periods r. Its variation with the 
control parameter 7 ex is shown in Fig. 0. The error bars 
represent the variance of the chaotic fraction. 
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FIG. 6: x-t-plots at different states of the data processing, (a) 
raw data; (b) same section as in (a) after transforming the 
data. White corresponds to large local wavelengths, black 
to small ones, (c) final step: the relative change of local 
wavelength over two subsequent periods A is calculated and 
binarized. The black areas are defined as regular (A < 0.01), 
the white are chaotic (A > 0.01). 
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FIG. 7: The mean chaotic fraction 7 vs. excitation amplitude 
/ex- The solid line is a power law fit. The error bars represent 
the statistical errors. 
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FIG. 8: The mean laminar expansion in space £ and in time 
T vs. control parameter e. The lines are power law fits. The 
error bars represent statistical errors. 



Close to the onset of STI the mean chaotic fraction is 
expected to grow with a power law 



7 ~ (te- 



le 



(2) 



The solid line in Fig.0 is a fit to our data, using 7 C , /?, 
and an offset representing background noise as adjustable 
parameters. The threshold value determined in this way 
is J c = 3.0 ± 0.05 A and the exponent /3 = 0.3 ± 0.05 
is in agreement with the theoretical expectation for DP 
(3 = 0.276486(8) |§. 

Another way to characterize the regular domains is to 
look at the mean laminar expansion in space £ and time 
T. First we define the laminar lengths I as the number of 
consecutive regular segments between two chaotic ones 
devided by the total number of segments at a certain 
time to- The laminar times t are the number of segments 
between two chaotic ones at a certain position x$. The 
averages of these numbers are displayed in Fig.|§] as the 
function of the normalized controlparameter 



= — 1 



Both parameters decay with a power law: 



T ~ e~ 



(3) 



(4) 
(5) 



and the fit yields the exponents v* = 1.2 ± 0.2, which is 
represented by the solid line and v* = 0.7 ± 0.1, which is 
represented by the dashed line. Only data in the range 
of 0.03 < e < 0.1 are taken into account. For smaller e a 
finite size effect is obvious: Following Cross and Hohen- 
berg |29| the characteristic length of the regular domains 
£ must be much smaller than the system size 1. For 
e > 0.1 the system is no longer intermittent, but rather 
chaotic. 
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FIG. 9: Distribution of the laminar domain length / and time 
t for e = 0. The solid lines are power law fits. To suppress 
the statistical fluctuations the distributions are logarithmical 
binned. 



In Fig. g the distributions of the laminar domain length 
and time for e = are presented. At the threshold the 
distribution should follow a power law for both the dis- 
tribution of the laminar domain lengths I and times t. To 
suppress the statistical fluctuations the values are loga- 
rithmically binned. The solid line represent the power 
law fit for the distribution of the lengths (||) and for the 
times ([?]): 



p(0< 

p(t) 



(6) 
(7) 



We obtain /i s = 1.7 ± 0.05, in agreement with the theo- 
retical value /i s = 1.734 and fi t = 2.1 ± 0.1 in accordance 
with the theoretical value /i t > 2.0. 

For e > the power law has a cutoff at the correlation 
length resp. time of the system and an exponential tail 
with a decay length equal to the correlation length £decay 
resp. time (9 deC ay (Fig.[l(]). 

Both parameters should grow with the same power law 
as the mean expansions £ and T in Eqs. (|4|) and (||). This 
behaviour can be seen in Fig.|ll|: The lines correspond 
to power law fits with the exponents u s = 1.1 ±0.2 and 
Vt = 0.62±0.14. The large errors are due to the statistical 
fluctuations of the distributions and difficulties in the 
definition of the cutoff length. 

The distributions cannot be described by a simple 
power law. A more complicated distribution function 
has been suggested 



p(l) = (M-* + B)e 



-l/lde- 



(8) 



in Ref. fllql . The solid line in Fig. [12] is a fit to this empiri- 
cal distribution function for e = 0.019 with /decay = 0.17. 
It shows clearly that the power law is now replaced by a 
function more reminiscent of an exponential decay. 
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FIG. 10: Distributions of the laminar length / (a) and the 
laminar time t (b) for e > 0. The legend in (b) holds for both 
figures. 
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FIG. 11: The correlation length £ and time # derived from the 
exponential decay of the distributions (Fig. vs. control 
parameter e. The lines are power law fits. The error bars 
represent statistical errors. 
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FIG. 12: Distribution of the laminar domain length / for e 
0.019. The solid line stems from Eq. (^). 



IV. DISCUSSION 

To conclude, we have presented a new experimental 
system exhibiting STL In contrast to all previous experi- 
ments, which are autonomous ones, our system is period- 
ically driven. We have measured the critical exponents 
/?, z/g, ^t, /i s and fit for the mean chaotic fraction 7, the 
mean laminar length £, the mean laminar time T, the 
correlation length £decay and time #decay and the laminar 
length and time distribution functions for e ~ 0. Four 
of the five exponents agree with the theoretical expecta- 
tion derived from a DP model within our experimental 
resolution. Considering the simplicity of the underly- 
ing discrete model, the fact that the theory is applicable 
only near 7 C , the complexity of our experiment and the 
fact that our apparatus has a finite size this concordance 
seems truly remarkable. 

The fact that the parameter v t is far from the expected 
value needs further discussion. One difference between 
the underlying model and our experiment might be the 
nature of the absorbing state. In fact, the absence of a 
truly absorbing state seems to be common to all experi- 
mental results so far [|o|. Following the ideas presented 
in that book the allowance for a continuum of states be- 
tween the absorbing and the active one might result in a 
more realistic description of the experimental situation. 

Another way to soften the assumption of a truly ab- 
sorbing state is the introduction of stochastic mecha- 
nisms permitting the nucleation of chaotic domains. In 
DP this is done by implementing a weak external field 
which creates chaotic domains in a pure laminar neigh- 

. As long 



bourhood with a certain probability [23, | 
as the probability of the creation processes is small the 
critical behaviour of DP is only slightly disturbed. If 
the probability becomes too large the universality is de- 
stroyed. The stochastic field could be an anologon to the 
"background noise" of our experimental apparatus. 
Alternatively, traveling soliton-like structures have 
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added to CML and DP [|2|. The interaction of these 
structures might lead to chaotic domains nucleating in 
laminar regions, which again softens the assumption of 
a truly absorbing state and leads to the breakdown of 
universality in DP. 

If any of these ideas will be able to explain the be- 
haviour of the measured decay times is subject to further 
investigation. 
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